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Abstract 



The adsorption of atomic oxygen and its inclusion into subsurface sites on Ag(210) 
and Ag(410) surfaces have been investigated using density functional theory. We 
find that — in the absence of adatoms on the first metal layer — subsurface adsorption 
results in strong lattice distortion which makes it energetically unfavoured. However 
subsurface sites are significantly stabilised when a sufficient amount of O adatoms is 
present on the surface. At high enough O coverage on the Ag(210) surface the mixed 
on-surface + subsurface O adsorption is energetically favoured with respect to the 
on-surface only adsorption. Instead, on the Ag(410) surface, at the coverage we 
have considered (3/8 ML), the existence of stable terrace sites makes the subsurface 
O incorporation less favourable. These findings are compatible with the results of 
recent HREEL experiments which have actually motivated this work. 

Key words: Density functional calculations. Catalysis, Silver, Oxygen, Stepped 
single crystal surfaces. 



1 Introduction 

The interaction between oxygen and silver surfaces has been the subject of ex- 
tensive experimental and theoretical research because of its key role in impor- 
tant heterogeneous catalytic reactions, such as partial oxidation of methanol to 

* Corresponding Author. Phone: +39-040 3787477, Fax: +39-040-3787528. 
Email address: bonini@sissa.it (N. Bonini^). 



Preprint submitted to Elsevier Science 



2 February 2008 



formaldehyde or ethylene epoxidation [1]. Despite significant efforts, these sys- 
tems are still not well understood and the identification of the active species 
involved in these reactions remains unclear. In particular, defect sites and 
subsurface oxygen species have attracted considerable attention since they are 
believed to be important in the ethylene epoxidation reaction [2]. 

In recent experimental works [3,4,5] Rocca and coworkers have investigated 
the interaction of oxygen with Ag(410) and Ag(210) surfaces with the aim 
of understanding the role of steps in the adsorption of oxygen on silver. In 
particular, using a supersonic molecular beam to dose oxygen on the sur- 
face and characterising the final adsorption state by vibrational spectroscopy 
(HREEL), they have studied the different oxygen species which form on these 
surfaces. Of particular interest are the results for O/Ag(210) [5]. The HREEL 
spectra of this system show a peak around 56 meV which is proposed to be 
due to the vibration of an oxygen atom occupying a subsurface site. Such a 
peak appears only when preadsorbed oxygen is present on the surface and it 
is not observed neither on low Miller index surfaces nor on O/Ag(410). 

Motivated by these results we have performed an ab initio calculation to inves- 
tigate the possibility of incorporating oxygen into interstitial subsurface sites 
on Ag(210) and Ag(410) surfaces. 

Recent theoretical works have studied the inclusion of O into subsurface re- 
gions of low Miller index silver surfaces, i.e. Ag(lll) [6] and Ag(lOO) [7]. It 
was found that pure subsurface O atoms are always unstable with respect 
to the oxygen molecule but the presence of co-adsorbed on-surface O atoms 
can stabilise them. More generally, we find that the adsorption in subsur- 
face sites depends significantly on coverage. Isolated subsurface oxygen atoms 
are always less stable than adatoms occupying favourable on-surface sites, 
essentially because of the high energy cost for lattice distortion induced by 
subsurface O incorporation. However, increasing the coverage the preference 
of on-surface sites decreases and, for example, in 0/Ag(lll) pure subsurface 
sites are stabilised at coverages higher than 0.5 ML [6]. This behaviour is due 
to the fact that the repulsive interaction between the adsorbates is more effec- 
tively shielded for the subsurface atoms than for the on-surface ones. In this 
work we consider high Miller index surfaces in order to investigate the effects 
of steps on the incorporation of oxygen atoms into subsurface sites. 



2 Computational Details 

Calculations were performed in the framework of density functional theory 
(DFT) using the generalised gradient approximation (GGA) of Perdew-Burke- 
Ernzerhof (PBE) [8] . We have used the pseudopotential method with ultra-soft 
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pseudopotentials [9] and plane- wave basis sets up to a kinetic-energy cutoff of 
27 Ry (216 Ry for the charge-density). Details about the Ag and O pseudopo- 
tentials are reported in Ref. [10]. Brillouin zone (BZ) integration has been 
performed with the Gaussian-smearing special-point technique [11,12] with a 
smearing parameter of 0.03 Ry. Calculations have been done using the PWscf 
package [13], while the molecular graphics were produced by the XCRYSDEN 
[14] graphical package. 

We have used periodic super-cells as models of the surface. The Ag(210) sur- 
face is modelled with slab of 14 (210) layers, which corresponds to seven (100) 
layers rotated by an angle 02 (tan(02) = 1/2) around the [001] axis and merged 
together so as to form the Ag(210) surface. The Ag(410) is modelled with slabs 
of 20 (410) layers corresponding to five (100) layers rotated by an angle 04 
(tan((/)4) = 1/4) and merged together. Adjacent slabs are separated by a vac- 
uum region of at least 16 a.u. Oxygen is adsorbed on both sides of the slab 
and all the structures have been fully relaxed until the Hellmann-Feynman 
forces were lower than 10~^ Ry/a.u. 

Adsorption of atomic oxygen has been modelled by (2x1) surface super-cells. 
For the Ag(210) surface the BZ integrations have been performed using a 
(4x7x1) Monkhorst-Pack mesh, while for the Ag(410) surface a (4 x 4 x 1) 
mesh has been applied. All the calculations are non-spin-polarised, because 
the interaction between O and the Ag substrate results in a negligible spin 
moment on the O atom. 

The chemisorption energies, E'chem, are referred to the clean Ag(nlO) sur- 
face (n = 2,4) and the isolated oxygen molecule, Echem = (£'0/Ag — £'Ag — 
No(Eo^/2))/No^ where the total energy of the adsorbate-substrate system, 
of the clean surface, of the isolated O2 molecule, and the number of adsorbed 
oxygen atoms are represented by ^^o/Ag, ^Ag, Eq^ and A^o^ respectively. £^02 
is estimated from a spin-polarised calculation. With this definition, stable ad- 
sorbates have negative chemisorption energies. 



3 Results and Discussion 

In previous works on Ag(OOl) [10,15], we found that the surface four-fold hol- 
low site is the most stable one and the chemisorption energy was found to 
decrease increasing the coverage, a fact related to the repulsive lateral elec- 
trostatic interaction between negatively charged oxygen adatoms. Recently, 
we performed a systematic study of on-surface atomic oxygen adsorption on 
Ag(410) and Ag(210) surfaces [16]. Figure 1 shows the structure of the two 
surfaces and the on-surface sites that have been studied. In site A the oxy- 
gen adatom lies just above the step and is coordinated with three Ag surface 
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atoms; in site B the oxygen adatom lies in the hollow site just below the step. 
These two sites are present on both surfaces. The Ag(410) surface possesses 
also the Ti and T2 sites (see Figure 1) with the O adatom in the four- fold 
hollow sites on the (100) terrace. We find that the configuration in which the 
adatoms occupy the sites A forming -0-Ag-O-Ag- rows (configuration la- 
belled A-A) is particularly stable. On the Ag(210) surface the step decoration 
is significantly more favourable than the other adsorption configurations. On 
Ag(410), instead, the terrace sites have a chemisorption energy almost de- 
generate with that of the step sites. In particular, on this surface, when the 
adatoms are far apart from each other the hollow sites on terrace are slightly 
more stable than the step edge sites, while at higher coverage the adatoms 
slightly prefer to decorate the steps. 




Fig. 1. Schematic representation of the Ag(210) (left) and Ag(410) surface (right). 
The labels A, B, Ti and T2 indicate different adsorption sites: A, above the step; 
B, below the step; Ti and T2, hollow sites on terrace. The surface unit-cell vectors, 
and the [001] and [InO], n = 4, 2, crystal axes are also shown. 

Let us focus now on the energetics of oxygen inclusion into subsurface intersti- 
tial sites, addressing in particular the effects of on-surface co-adsorption. Our 
results are summarised in Table 1 and discussed in the next paragraphs. 



The inclusion of oxygen into subsurface sites is strongly unfavourable with 
respect to on-surface O adsorption, the resulting chemisorption energies on 
Ag(210) surface being as low as —0.07 eV for subsurface-octahedral site (Octa) 
and —0.14 eV for subsurface-tetrahedral site (Tetra).^ It is interesting to 
observe that the adsorption in OCTA and Tetra sites on Ag(lOO) at ^ = 
0.25 ML [7] is considerably less stable (0.79 eV and 0.69 eV, respectively) 
than on Ag(210) surface. This indicates that the presence of steps induces a 
stabilisation of subsurface oxygens. In the right panels of Figure 2 we show 
top views of the optimised structures with O atom in the OCTA and Tetra 
sites. Note that due to a significant distortion of the substrate lattice in- 
duced by subsurface O atoms, the surface structure of Ag(210) can hardly be 
recognised. However these subsurface sites are significantly stabilised when a 

^ These calculations were performed with Ag slabs composed of only 10 (210)-layers 
and somewhat smaller (3x5x1) Monkhorst-Pack mesh. 
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Table 1 

Average chemisorption energies, Echem^ for purely subsurface, mixed on-surface 
+ subsurface, and on-surface only adsorption of atomic oxygen on Ag(210) and 
Ag(410) surfaces. The label of OCTA, Tetra, A, B, and Ti are defined in text. 



Surface 


Coverage 


Configuration 


^chem(eV) 


Ag(210) 


e = 1/4 ML 


OCTA 

Tetra 
A 


-0.07 
-0.14 
-0.68 




e = 3/4 ML 


A-A + OCTA 
A-A + B 


-0.54 
-0.40 


Ag(410) 


e = 3/8 ML 


A-A + OCTA 

A-A + B 
A-A + Ti 


-0.55 
-0.68 
-0.63 



sufficient amount of O adatoms is present on the surface. The chemisorption 
energies for A-A+OCTA configuration^ are —0.54 and —0.55 eV for Ag(210) 
and Ag(410) surfaces, respectively. Also observe that the substrate reconstruc- 
tion is considerably reduced in these cases, and that the Ag(210) and Ag(410) 
surface structures can be easily recognised from the corresponding optimised 
0/Ag structures (see left panel of Fig 3 and right panel of Fig 4). 

It is interesting to compare the stability of mixed on-surface + subsurface con- 
figurations with the purely on-surface ones at a given (high enough) oxygen 
coverage. In all the configurations we have considered, O adatoms decorate the 
step-edge (A-A arrangement). An additional oxygen atom is then considered 
to be located either into subsurface OCTA site (A-A+Octa) or adsorbed 
on a given surface hollow site (see Figs 3 and 4). For Ag(210), we consider 
only the A-A-|-B configuration with an additional O adatom located in the 
on-surface hollow B site situated below the step edge. For Ag(410), in ad- 
dition to the A— A+B geometry, we also consider the configuration with an 
additional O atom adsorbed on the terrace Ti site (configuration A-A+Ti). 
On Ag(210) the mixed on-surface + subsurface A-A+OCTA configuration 
is more stable than the pure on-surface A-A+B configuration, Echem being 
—0.54 and —0.40 eV, respectively. On the contrary, the opposite is true for the 
Ag(410) surface, where the chemisorption energies for the A— A+OCTA and 
A-A+B configurations are —0.55 and —0.68 eV, respectively. This qualita- 
tive difference is related to the different O coverages of the two surfaces, which 
is larger on Ag(210) surface {9 = 3/4 ML) than on Ag(410) {9 = 3/8 ML). 



^ In this configuration O adatoms decorate the step-edge and an additional O atom 
per (2 X 1) super-cell is located into OCTA site situated straightly below a step-edge 
Ag atom. 
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atom in octahedral- subsurface site 

Initial structure Optimised structure 




atom in tetrahedral- sub surface site 

Initial structure Optimised structure 




Fig. 2. Oxygen atoms located in subsurface sites on Ag(210) surface. Large (small) 
balls are silver (oxygen) atoms. On left-side panels the horizontal positions of subsur- 
face-octahedral (Octa) (top left panel) and subsurface-tetrahedral (Tetra) (bot- 
tom left panel) sites are marked with crosses (x), while vertically these sites are 
located below the first layer Ag atoms in the octahedral and tetrahedral interstices. 
On right-side panels the corresponding optimised 0/Ag structures are shown. Note a 
significant distortion of the substrate lattice — the (210) surface-structure attributes 
can hardly be recognised. 



A-A ^i^A 





Fig. 3. The two competing O/Ag(210) structures at ^ = 3/4 ML oxygen coverage. 
On both structures the oxygen atoms are decorating the step-edge (A-A). On the 
left structure an additional oxygen per (2 x 1) super-cell is located in the subsur- 
face-octahedral site (Octa), i.e., O atoms located straightly below the marked (x) 
step Ag atoms in octahedral interstices. On the right structure the additional O 
atom per (2 x 1) super-cell is located in the hollow site below the step-edge (B). 
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Fig. 4. The O/Ag(410) structures at ^ = 3/8 ML oxygen coverage. From left to 
right: A-A + OCTA, A-A + B and A-A + Ti. For the A-A + OCTA structure 
the OCTA-oxygen atoms are located straightly below the marked (x) step Ag atoms 
in octahedral interstices. 

This suggests that the electrostatic repulsion between the adsorbates in the 
A-A+B configuration is stronger on Ag(210) than on Ag(410) and makes 
this geometry on Ag(210) significantly less favourable. Also note that the 
eflFect of high coverage in Ag(210) is weaker for the A— A+OCTA configura- 
tion essentially because the distance between the adatoms is larger than in 
the A— A+B arrangement and because the interaction involving subsurface 
O species is shielded more effectively than for purely on-surface O atoms. In 
addition to this, we have found that on Ag(410) the incorporation of oxygen 
into subsurface sites is less stable than on-surface adsorption also because of 
the availability of stable terrace sites. The A— A+Ti configuration is indeed 
more stable than the A— A+OCTA one. 



4 Conclusions 

Our results are compatible with the suggestion, put forward in Ref. [5], that 
the HREELS peak found at 56 meV in Ag(210), and missing in Ag(410), is 
due to atomic oxygen adsorbed subsurface. Indeed our calculations, although 
performed at different coverages, indicate that subsurface O is not stable in 
Ag(410), while it is stabilised in Ag(210) at high enough coverage. Further 
theoretical investigations spanning a larger range of coverages are called for. 
Lattice vibrational calculations on selected 0/Ag(nlO) configurations would 
probably fully clarify the picture. 
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